MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BURIAU  OF  ST ANDARDS- 1963  A 


DOC  FILE  COPY  b ADA037  3 72 


MOST  Project  - -z- 


v 


THE  UNIVERSITY  OF  TEXAS 

AT  AUSTIN 


July  1971  c°py  No-  ' 

QUARTERLY  PROGRESS  REPORT  NO.  1 

UNDER  CONTRACT  N00024-7 1-C- 1 266  - 

April,  May,  and  June  1971 


NAVAL  SHIP  SYSTEMS  COMMAND 

Contract  N0"024-71  -C-l  266 

Proj.  Ser.  N0.  SF  11552002,  Task  8118 


D D C 


HrG07-CC75 


0 * 
SQt 


DISTRIBUTION  STATEMENT-  A 

Approved  for  public  release 
Distribution  Unlimited 


Ksm  i«> 


I WSIk  SKtiN 

tor.  futin  o 

V.V.u^eSi  q 

. I 


•V  B')P£S 


fi 


22  July  1971 


QUARTERLY  PROGRESS  REPORT  HO.  1 
UNDER  CONTRACT  N00024-71-C-1266 

April,  May,  and  June  1971 


NAVAL  SHIP  SYSTEMS  COMMAND 

Contract  N00024-71 -C-l  266 

Proj.  Ser.  N0.  SF  11552002,  Task  8118 


D D C 


APPLIED  RESEARCH  LABORATORIES 

THE  UNIVERSITY  OF  TEXAS  AT  AUSTIN 

AUSTIN,  THXAS  78712 

DISTRIBUTION  SffiTEMgNT-  ' 

Approved  fai  public  release) 

Distribution  Unlimited 


ABSTRACT 


In  Section  I a status  report  is  given  on  the  experimental  measurements 
of  forward  scattering  by  acoustically  penetrable  rough  surfaces. 

Plastic  molds  from  the  existing  pressure  release  surfaces  have  been 
made.  One  penetrable  rough  surface  has  been  cast.  The  other  surfaces 
will  be  completed  as  soon  as  the  acoustical  parameters  of  the  first 
one  have  been  studied.  In  Section  II  the  acoustic  waveguide  problem 
for  lossy  boundaries  is  solved  in  terms  of  a contour  integral.  The 
eigenvalue  equation  is  given  as  a function  of  the  reflection  coefficients 
of  the  lossy  boundaries. 
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I.  PRELIMINARY  WORK  ON  THE  EXPERIMENTAL  MEASUREMENTS  OF 
FORWARD  SCATTERING  BY  ACOUSTICALLY  PENETRABLE  SURFACES 

In  preparation  for  the  experimental  forward  scattering  measurements 
of  acoustically  penetrable  random  surfaces,  the  existing  pressure  release 
surfaces  were  used  to  make  plaster  molds.  From  these  molds  surfaces 
of  the  same  statistical  properties  as  the  pressure  release  surfaces  but 
different  acoustical  properties  are  to  be  cast.  Initial  plans  called 
for  molds  of  the  smoothest  and  the  roughest  of  the  three  surfaces 
(rms  heights  = 0.091  in.,  0.182  in.,  and  0.36U  in.  and  identical  cor- 
relation lengths)  to  be  made.  However,  the  roughest  surface  could 
not  be  used  for  mold  making  because  of  deterioration  of  the  surface 
coating.  Consequently,  plaster  molds  of  only  the  smoothest  and  inter- 
mediate surface  were  made.  Casting  of  these  surfaces  using  an 
acoustically  penetrable  material  is  in  progress.  In  addition,  an 
acoustically  penetrable  plane  surface  will  be  made  for  data  comparison. 

The  material  which  will  be  used  to  simulate  a liquid  bottom  (low 
shear  wave  velocity  bottom)  is  a mixture  of  Scotchcast  221  polyurethane 
and  No.  5 sandblasting  sand  in  the  proportion  1:1  by  weight.  The 
selection  of  suitable  materials  and  the  measurement  of  their  acoustical 
properties  were  carried  out  during  the  previous  contract  year  and  are 
reported  in  detail  in  Appendix  A of  the  Final  Report  under  Contract 
N0002U-70-C-1279*  "Scattering  and  Propagation  of  Acoustic  Waves  in  the 
Presence  of  Rough,  Penetrable  Boundaries."  A penetrable  surface  has 
been  cast  from  the  mold  of  the  smoothest  surface  (rms  height  = 0.091  in.) 
however,  further  casting  has  been  temporarily  halted  to  allow  for  an 
evaluation  of  the  acoustic  parameters  of  the  new  surface.  This  is 
necessary  because  of  the  unforeseen  settling  of  the  sand  in  the  casting 


material.  The  original  check  of  the  acoustical  properties  used  small 
volumes  of  material  and  showed  no  evidence  of  settling.  However, 
during  the  recent  casting,  large  volumes  of  material  were  used  and 
settling  of  the  sand  was  observed.  Since  settling  of  the  sand  creates 
a density  gradient  in  the  surface,  a casting  method  which  maintains 
the  homogeneity  of  the  mixture  is  presently  being  sought. 

Upon  completion  of  the  surface  castings  and  evaluation  of  the 
acoustic  parameters  of  these  surfaces,  the  experimental  forward 
scattering  measurements  will  begin. 


II.  NORMAL  MODE  SOUND  PROPAGATION  IN  AN 
INHOMOGENEOUS  MEDIUM  WITH  LOSSY  BOUNDARIES 


A.  Introduction 

The  acoustic  field  of  a cv  point  source  located  in  an  inhomogeneous 
medium  with  lossy  boundaries  is  given  as  a contour  integral.  Singu- 
larities of  the  contour  integral  are  given  by  the  dispersion  equation 
(period  equation,  eigenvalue  equation,  etc.).  Lossy  boundaries  are 
introduced  through  effective  reflection  coefficients.  By  applying 
phase  integral  techniques,  the  dispersion  equation  is  given  in  terms 
of  the  reflection  coefficients  of  the  lossy  boundaries. 

B.  Fundamental  Integral  for  the  Acoustic  Field 

The  wave  equation  for  a point  source,  with  angular  frequency  w, 
located  in  a layered  inhomogeneous  medium  is  given  by 
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where  P represents  the  pressure,  C(z)  is  the  sound  velocity  (variable 

in  the  depth  coordinate  z) , and  6(r-?Q)  is  the  three-dimensional  Dirac 

delta  function.  The  strength  Q of  the  source  is  equal  to  the  density 

o 

times  the  volume  velocity.  If  cylindrical  coordinates  (r,9,z)  are 
assumed  with  azimuthal  symmetry  and  a time  factor  of  exp(-iut)  is 
suppressed,  the  following  equation  is  obtained: 
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where  P now  represents  the  time  independent  pressure.  The  wave  number 

k(z)  is  defined  as  w/C(z)  and  the  source  is  located  at  z=z  and  r® 0, 

o 

as  indicated  in  Dwg.  AS-70-767.  The  vertical  depth  coordinate  z varies 
from  0Sz<d,  and  the  range  coordinate  r varies  from  0<r<°°.  Boundary 
conditions  imposed  on  Eq.  (2)  are  that 


a)  P must  satisfy  a radiation  condition  for  r-*» 

b)  at  the  surface  and  bottom,  the  reflection  coefficients  be  known 

3P 

c)  P and  — be  continuous  across  any  discontinuities  in  the 

3 z 

velocity  depth  profile. 

By  means  of  a convolution  of  Green's  functions  the  solution  for 
Eq.  (2)  can  be  stated  as 
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where  £ is  the  complex  separation  parameter,  and  where  z<  and  z>  denote 
the  smaller  or  larger,  respectively,  of  the  variables  z and  zq.  The 
functions  ^ and  \\>2  are  chosen  to  satisfy  boundary  conditions  at  z=0 
and  z=d.  When  boundary  conditions  are  other  than  pressure  release  and 
rigid,  the  usual  method  is  to  state  the  boundary  values  as  impedance 
conditions.  For  example  Bucker1  (in  his  Appendix  A)  states  the  boundary 
conditions  as  impedance  conditions.  By  means  of  a pseudoisovelocity 
layer,  Bucker  relates  his  impedance  conditions  to  plane  wave  reflection 
coefficients.  He  assumes  various  values  for  the  boundary  losses  and 
calculates  the  reflection  coefficients  from  the  equations 
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where  SL  and  BL  represent  surface  loss  and  bottom  loss  respectively, 
and  where  S and  R are  the  surface  and  bottom  plane  wave  reflection 
coefficients.  Therefore  by  assuming  various  values  for  the  surface 
and  bottom  losses,  Bucker  is  able  to  arrive  at  Impedance  conditions 
at  the  boundaries. 

C.  The  Wronskian  Dispersion  Equation  • 

According  to  Appendix  A of  Bucker 's  paper  the  dispersion  equation 
is  given  by  setting  the  Wronskian  equal  to  zero,  that  is, 

W(5  ) = 0 (5) 

n 

However,  Eq.  (5)  does  not  always  give  the  correct  results.  For 
example,  in  the  case  of  surface  duct  propagation  where  no  bottom  is 
assumed  (i.e.,  z-h»)  and  where  an  impedance  condition  is  assumed  at 
the  surface,  then 


^(z)  = n1(z,C)  -X  n2(z,C) 

and 


(6) 


^(z)  = n2(z,C)  , (7) 

where  n^  and  n^  are  traveling  wave  solutions  to  the  z-component  wave 
equati on : 


(8) 


Solution  n?  must  represent  an  outgoing  wave  for  z-**>  and  n^  represents 
an  outgoing  wave  for  z*-®. 
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The  function  X is  given  by 


where  the  impedance  condition  is  given  by 


In  this  case  the  eigenvalues  should  be  given  by  the  zeros  of  the 
Wronskian  of  i|>  and  : 


but  the  correct  period  equation  is  given  by 


It  appears  that  no  general  form  for  the  secular  equation  can  be  given 
in  terms  of  the  Wronskian.- 


D.  The  Reflection  Coefficient  Dispersion  Equation 


The  eigenvalue  equation  in  terms  of  effective  reflection  coefficients 
is  given  by  Brekhovskikh,2  Budden,3  and  Wait.4  The  general  form  is 
stated  as 


2ikh  sin  0 


where  the  region  between  zero  and  h represents  a homogeneous  layer. 

The  regions  above  z=0  and  below  z=h  may  be  inhomogeneous.  and 
are  each  functions  of  the  modal  parameter  £.  A form  for  (bottom 

reflection)  has  been  given  by  Urick,  Lund,  and  Bradley.5  From  classical 
scattering  theory  the  form  for  V (surface  reflection)  can  be  given  as 


v1(e)  = -e 


2(ko  sin  0)‘ 


where  o is  the  rms  height  value  for  the  rough  ocean  surface.  If  a 
peak  to  trough  height  H is  known,  then  a good  approximation  for  a is 
given  by  the  equation 


(15) 


It  is  interesting  that  no  suitable  forms  for  an  eigenvalue  dependent 
impedance  condition  are  found  in  the  literature,  a suitable  form  being 
one  which  is  both  surface  property  dependent  and  medium  dependent. 
However,  if  the  boundary  conditions  are  specified  in  terms  of  reflection 
coefficients,  then  both  the  surface  properties  and  the  medium  properties 
(through  the  angle  fl)  are  included. 
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For  the  situation  where  there  is  no  homogeneous  region,  the 
parameter  h is  allowed  to  approach  zero.  V and  are  then  effective 
upgoing  and  downgoing  reflection  coefficients  referred  to  the  same 
level.  Since  the  Wronskian  evaluation  does  not  depend  on  any  specific 
level,  neither  does  the  level  for  and  V . Equation  ( 13 ) then  becomes 

v V = l (16) 

So  long  as  there  existed  a homogeneous  layer  between  zero  and  h , 
the  lossy  boundary  value  problem  in  terms  of  effective  reflection 
coefficients  was  straightforward.  When  the  medium  is  inhomogeneous, 
then  it  is  not  so  clear  how  to  include  both  boundaries  when  the  two 
reflection  coefficients  are  referred  to  the  same  level.  However  if 
phase  integral  methods  are  used  to  obtain  the  period  equation,  then 
the  reflection  coefficients  can  be  referred  to  the  two  boundaries. 

The  phase  integral  period  equation  should  give  fair  estimates  for  the 
generally  complex  eigenvalues.  Bartberger  and  Ackler8  have  discussed 
the  value  of  the  phase  integral  condition  for  the  Furry-Pedersen  type 
of  velocity  depth  profiles. 

The  phase  integral  period  equation  can  be  stated  as7’8 


V1(o)  V?(d) 


1 


(IT) 


which  1 s a generalization  of  the  usual  condition  without  the  V's. 

It  is  instructive  to  write  Eq.  (IT)  in  a transverse  resonance  form: 


dz  = Pmr  + iHn  V (oJV^fd) 


(n*0 ,1,2,...) 


(18) 
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This  form  is  essentially  the  same  as  Bucker9  derived  several  years  ago. 
The  usual  phase  integral  conditions  can  be  obtained  from  Eq.  (18). 

For  example  if  V^o)  = e“l7r  = -1  (free  surface)  and  V^fd)  = e 11 
(mode  trapped  by  continuous  refraction  due  to  the  inhomogeneous  medium 
at  some  point  d) , then  Eq.  (l8)  becomes 


dz  = ir(n  + 3/1*) 


(n=0,l,2 , . . . ) 


(19) 


When  the  modes  are  completely  trapped  because  of  continuous  refraction 

— i /2 

due  to  the  inhomogeneous  medium,  then  V^o)  = V^d)  = e w and 
Eq.  (18)  becomes 


J ^ k2(z)  - £2  dz  = n(n  + 1/2)  . (20) 

It  should  be  noted  that  when  the  modes  are  trapped  because  of  the 
medium,  then  the  lines  z*o  and  z=d  Eire  caustics  of  the  guided  ray 
system  and  the  n/2  values  give  the  phs^e  change  of  a ray  upon  touching 
the  caustic. 

By  separating  Eq.  (l8)  into  reed.  Euid  imaginary  parts  it  is  possible 

to  obtadn  forms  for  Re(r  ) and  Im( f ).  The  form  for  Im( £ ),  which  is 

^n  sn  n 

the  result  of  boundeiry  losses  Euid  leakage,  can  sometimes  be  expressed 
in  terms  of  equivalent  ray  cycle  distsuice.  These  values  for  ^ are 
only  approximate  and  will  only  be  used  for  estimates  to  the  exact  values. 
The  exact  values  must  be  obtained  from  Present  work  is  being 

directed  in  obtaining  an  alternate  form  for  * 1. 
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